What is claimed is: 

1. A numerical analysis system performing numerical analysis by Finite 
Volume Method employing Hybrid Grid Adaptation Method comprising: 

(a) an initial grid memory containing the initial grid data made by way of 
free combination of tetrahedral cells, hexahedral cells and pentahedral 
cells; 

(b) a grid dividing means dividing a cell into plural cells in order to divide 
a triangular face of a cell into plural triangular faces and divide a 
quadrilateral face of a cell into plural quadrilateral faces; and 

(c) a grid deleting means recovering the former cell by deleting divided 
cells. 

2. The numerical analysis system of claim 1, wherein said grid deleting 
means comprising: 

(cl) a dividing threshold memory containing the dividing threshold; 
(c2) a deleting threshold memory containing the deleting threshold; 
(c3) a computing means to compute the divide-delete index depending 
upon the density in each cell; 

(c4) a division deciding means deciding to divide the cell when the volume 
of the cell is greater than the minimum and when the division level of the 
cell is not maximal and when the divide-delete index is greater than the 
dividing threshold; and 

(c5) a deletion deciding means deciding to delete the cell when the cell is 
not the initial cell of zero division level and when the divide-delete index is 
less than the deleting threshold. 

3. The numerical analysis system of claim 1, wherein said grid dividing 
means comprising: 

(bl) a prism dividing means dividing a prismatic cell into 8 prismatic cells; 
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(b2) a hexahedron dividing means dividing a hexahedral cell into 8 
hexahedral cells; 

(b3) a pyramid dividing means dividing a pyramidal cell into 6 pyramidal 
cells and 4 tetrahedral cells; 

(b4) a simple tetrahedron dividing means dividing a tetrahedral cell into 8 
tetrahedral cells; and 

(b5) a tetrahedron dividing means dividing a tetrahedral cell into 4 
tetrahedral cells and 2 pyramidal cells. 

4. A numerical analysis method processing numerical analysis by Finite 
Volume Method employing Hybrid Grid Adaptation Method comprising: 

(a) an initial FVM processing step to calculate the quantity including 
density of fluid by Finite Volume Method depending upon the initial grid 
data made by using tetrahedral cells, hexahedral cells and pentahedral cells; 

(b) an index processing step to calculate the divide-delete index of each cell 
depending upon the calculated density; 

(c) a division processing step to divide a cell into plural cells in order to 
divide a triangular face into plural triangular faces and to divide a 
tetrahedral face into plural tetrahedral faces when the volume of the cell is 
greater than minimum and when the division level of the cell is not 
maximum and when the divide-delete index is greater than the division 
threshold; 

(d) a deletion processing step to recover the undivided cell by deleting such 
cells that the cell is not the initial cell of zero dividing level and the 
divide-delete index is less than the deleting threshold; 

(e) a renewal processing step to create the next grid; 

(f) an FVM processing step to calculate the physical quantities including 
the density of fluid by Finite Volume Method depending upon those grid 
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data; and 

(g) a repeating step to repeat those steps of (b) to (f). 

5. The numerical analysis method of claim 4, wherein a prismatic cell is 
divided into 8 prismatic cells, a hexahedral cell is divided into 8 hexahedral 
cells, a tetrahedral cell is divided into 6 pyramidal cells and 4 tetrahedral 
cells, a tetrahedral cell is divided into 8 tetrahedral cells or into 4 
tetrahedral cells and 2 pyramidal cells. 
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